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Abstract 

About 4 — 5% of chondrules are compound: two separate chondrules stuck together. This is commonly believed to be 
the result of the two component chondrules having collided shortly after forming, while still molten. This allows high 
velocity impacts to result in sticking. However, at T ~ HOOK, the temperature below which chondrules collide as solids 
(and hence usually bounce), coalescence times for droplets of appropriate composition are measured in tens of seconds. 
Even at 1025 K, at which temperature theory predicts that the chondrules must have collided extremely slowly to have 
stuck together, the coalescence time scale is still less than an hour. These coalescence time scales are too short for the 
collision of molten chondrules to explain the observed frequency of compound chondrules. We suggest instead a scenario 
where chondrules stuck together in slow collisions while fully solid; and the resulting chondrule pair was subsequently 
briefly heated to a temperature in the range of 900 — 1025 K. In that temperature window the coalescence time is finite 
but long, covering a span of hours to a decade. This is particularly interesting because those temperatures are precisely 
the critical window for thermally ionized MRI activity, so compound chondrules provide a possible probe into that vital 
regime. 
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1. Introduction 

Chondrules are sub-millimeter sized igneous inclusions 
found in chondritic meteorites which were melted in the 
Solar Nebula, the gas and dust disk from which our So¬ 
lar System formed. Usually distinct objects within their 
host meteorite, they are sometimes found as part of a com¬ 
pound chondrule: two distinguishable chondrules fused to¬ 
gether (see Figure [^. While compound chondrules (ccs, 
lower-case to distinguish from carbonaceous chondrites) 
are rare, with a frequency of about 4 — 5% of chondrules 


(Gooding & Keil, 1981 Ciesla et al., 2004), they are nonethe¬ 


less common enough for some basic statistical information 


about their nature to be well established (Gooding & Keil 


1981 Wasson et al. 19951. Wasson et al. (1995) distin¬ 


guishes between sibling and independent ccs, depending 
on whether the primary and secondary chondrules either 
have or don’t have similar composition and textural types. 
One of the greatest difficulties in studying compound chon¬ 
drules is the difficulty in correcting 2D information from 
thin sections to 3D properties (Ciesla et al. 2004), which 


makes determining their precise frequency difficult, and 
establishing other parameters, such as the size ratio of the 
secondary chondrule to the primary, fraught. 


Gooding & Keil (1981) proposed that compound chon¬ 


drules were made by collisions between still molten chon¬ 
drules shortly after they were formed from ambient dust. 
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If SO, then the cc frequency is an important constraint be¬ 
cause it links the chondrule-chondrule relative velocity and 


cooling times to the chondrule number density (e.g. Desch 


et al. 2012). However, sibling and independent ccs occur 


at similar frequencies, implying that if the cc fusing pro¬ 
cess occurred at temperatures where the chondrules were 
molten, that process must have been able to maintain the 
structural integrity of droplets of nearly identical liquids 
in close contact with each other. This requirement is ex¬ 
acerbated by the observation that small contact angle ccs 


dominate the statistics (Wasson et al. 1995 Ciesla et al. 


2004), i.e. the radius of the neck between the two compo¬ 


nent chondrules is generally much smaller than the radius 
of the individual chondrules, see Figure 


Wasson et al. (1995) suggests that at least some com¬ 


pound chondrules could have been made by melting neb¬ 
ular dust that had accreted onto a chondrule, which is a 
good model for enveloping ccs, a rare class of ccs so named 
because the second chondrule mostly envelopes the first. 


Miura et al. (2008) instead suggests a model where very 


large dust grains (far larger than observed chondrules) had 
their surfaces melted and stripped by a shock in the nebu¬ 
lar gas. The stripped surface separated into droplets which 
in turn collided with each other, creating compound chon¬ 
drules. 

Excepting the production of enveloping ccs by melting 
an accretionary dust rim, the above models all require that 
molten droplets remain in contact, with very small contact 
angles, for significant collisional time scales. Note that we 
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Table 1: H3 Chondrule Compositions 



BO 

Average 

SiOa 

45.4 

51.8 

TiOa 

0.21 

0.14 

AI 2 O 3 

5.7 

3.4 

Cr 203 

0.41 

0.55 

FeO 

14.0 

9.1 

MnO 

0.41 

0.34 

MgO 

28.7 

30.7 

CaO 

2.24 

2.30 

Na 20 

2.43 

1.40 

K 2 O 

0.45 

0.17 


Figure 1: A plane polarized light micrograph of a compound chon¬ 
drule from the Semarkona thin section AMNH 4128-1. 



Figure 2: Cartoon of a compound chondrule. The neck of the bridge 
has radius r, while (f> is the contact angle. 


use the term “molten” to refer to chondrules sufficiently 
liquid that they collide as liquids, rather than solids. This 
is important because it allows the individual chondrules 
to collide at sufficiently elevated speeds that the observed 
cc frequency might be matched. In this paper, we con¬ 
sider only non-enveloping ccs and show that even at tem¬ 
peratures well below chondrule formation temperatures 


(T > 1700K, Hewins & Radomsky, 1990), droplet coa¬ 
lescence time scales are short compared to collisional time 
scales and that the molten-collision model cannot match 
the observed cc frequency. That means that the small 
contact angles observed are a major constraint, and imply 
that cc fusing had to occur relatively cold. 


2. Viscosity and surface tension 

Chondrules exhibit a range of sizes, compositions and 


et al., 2004 


et al. 2014). Following previous studies of ccs (Ciesla 


textures (Weisberg 1987 Wasson et al. 1995 Friedrich 


Desch et al. 2012), we adopt a characteristic 


chondrule radius of a = 0.03 cm, and a chondrule solid 
density p = 3gcm“^. 

To understand the time scale on which molten chon¬ 
drules in contact with each other would have flowed we 


need their surface tensions and viscosities. Chondrule melts 
had surface tensions on the order of 


7 = 400 dyn cm 


( 1 ) 


or 0.4Nm Susa & Nakamoto 2002). 


We use Giordano et al. ( 2008 ) to calculate viscosities of 
chondrule melt^ Weisberg (19871 lists the composition of 
chondrules found in ordinary chondrites (see Table [^. We 
find that the viscosities for barred olivine/average chon¬ 
drules in ordinary chondrites were bracketed by the values 
of those chondrules in H3 chondrites: 


log ?7 = —3.55 -|- 


4557.8 
T-618.2 


and 


logr; = -3.55-h 


5084.9 


( 2 ) 


( 3 ) 


T - 584.9 

for barred olivines and average chondrules, respectively. 
Barred olivines were less viscous for relevant temperature 
ranges (see Figure left axis). T is the temperature 
in Kelvin and we are measuring the viscosity rj in poise 
(1 P = 0.1 Pa-s). While barred olivines are rare in the over¬ 
all chondrule population (Gooding & Keil 1981), using 


barred olivine compositions is not as strange as it might 
seem: Wasson et al. (1995) found that barred olivines 


are overrepresented in compound chondrules by a factor 
of many. 

There is clear evidence that chondrules experienced 


multiple heating events ( 

Jones & Danielson 1997 

Was- 

son & Rubin 

2003 

Jones et al. 

2005 

), and the existence 


of relict grains implies that at the temperatures we con¬ 
sider, many chondrules would have had more liquid, lower 
viscosity regions mixed with more solid, higher viscosity 


ones. However, as noted by Ciesla (2006), it is the lowest 


viscosity that controls the viscous behavior of the grains 
with the fluids flowing around the more solid inclusions. 


Ciesla et al. (2004) found that chondrules behaved as 


liquids while colliding if their interaction time was longer 


^At the time of writing, a convenient on-line cal¬ 
culator can be found at http://www.eos.ubc.ca/ krus- 
sell/VISCOSITY/grdViscosity.html 














































































































than the Maxwell time 


tm = (4) 

where /i ~ 4 x 10^° Pa is the shear modulus. They esti¬ 
mated the interaction time using Herztian contact theory, 
finding 


t,oii = 2.87 ( 

\mi+m2/ \ El E 2 J \ 0102 J 

(5) 

where m ~ 3.4 x 10 ^g, z/ ~ 0.25, and E ~ 10^^ Pa are 
the chondrules’ masses, Poisson ratios and Young’s moduli 
respectively. For those values, and a collision velocity Vc = 
100 cms“^, we can calculate 


tcoii ~ 3 X 10 ® s, ( 6 ) 

so the condition for liquid collision reduces to 

V<r]bc = 10® P. (7) 

That viscosity constraint implies 

T> 1095 77, ( 8 ) 

which value approximately applies for all collisional ve- 
locites as ry is a very steep function of temperature and 
the critical rjbc value is a very weak function of Vc- Below 
that viscosity limit, or equivalently above that tempera¬ 
ture limit, molten chondrules have a high sticking prob¬ 
ability even at velocities above 100cms“^. Below that 
temperature, chondrules collided as solids, bouncing unless 
the collisional speed was below O.lcms 

20iB. 


(Guttler et al. 


3. Coalescence time scales 


Liquid droplets in contact with each other tend to co¬ 
alesce into a single droplet. This makes it problematic 
that liquid-sticking requires 77 < r]bc = 10® P: we need the 
molten chondrules to have behaved as liquids for collisional 
interaction times, but not have behaved as liquids for far 
longer collisional or cooling windows. 

As might be expected for a phenomenon with signif¬ 
icant consequences for industry, there is a rich literature 


on droplet coalescence spreading ( 

Erenkel 

1945 

Tanner 

1979 

Eggers et al.l 1999[ |Yao et al. 2005 

Aarts et al. 


2005 


Eddi et al. 20131. Viscous droplet coalescence or 


spreading time scales depend on a coalescence time scale: 


arj 



(9) 


where a is the droplet radius (see Figure right axis). 
This scaling can be understood on dimensional grounds: 
larger and more viscous droplets take longer to coalesce, 
while the surface tension provides the energy source, and 



T[K] 


Figure 3: Viscosities (left axis) and coalescence times (right axis) for 
droplets with radii a = 0.03 cm. Solid/blue: H3 barred olivine chon- 
drule composition. Dashed/green: H3 average chondrule composi¬ 
tion. Vertical lines denote the temperature below which colliding BO 
chondrules bounce, and the temperatures at which BO chondrules 
have coalescence times r of ten seconds, an hour, a day and a week. 
The shaded area separates the viscosities below which droplets are 
sufficiently molten to stick well from the viscosities above which the 
individual droplets are sufficiently solid to maintain their structural 
integrity for meaningful time scales. 


T is the combination of a, r] and 7 with the correct di- 
mensions. One can also compare the energy dissipated by 
Stokes drag Fs on a sphere of radius a moving a distance a 
at a speed v through a medium of viscosity rj to the surface 
energy of a droplet with surface area ~ 47 ra^ and surface 
tension 7 : 

aFs = a X Gnriav ^ dTrya^, (10) 

and then estimate 


a arj 
r ^ ~ —. 

V 7 


( 11 ) 


Calculating r for 77 = rjbc and a = 0.03 cm, we find t = 75 s. 

The precise time dependence of the contact angle (or 
neck width) between two droplets remains a topic of re¬ 
search, but it is clear that there are multiple regimes where 
the radius of the interface between the droplets depends 
on t" for varying n, possibly with logarithmic corrections. 
However, it is nonetheless also clear from the literature 
that small contact angles can only be preserved for time 
scales 

t < r, ( 12 ) 


which places an upper limit on how long molten chon¬ 
drules can collide for. While this condition has the op¬ 
posite implication as the plasticity time tpias of previous 
work (Gooding & Keil 1981 Ciesla et al. 20041, it enters 
into analyses of cc collision rates in the same fashion. 

We should note that surface diffusion and sintering 
provides an alternate route to coalesence, but one with 
a characteristic time that scales as a^, significantly more 
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strongly than fluid coalesence with r oc a (Eggers 19981. 


As such, it is expected to act relatively slowly on the scale 
of a ~ 0.03 cm chondrules and indeed, extrapolating from 


Poppe (2003)’s data for Si 02 spheres, surface diffusion acts 


far slower than liquid coalescence for the chondrules we are 
interested in. 


4. Collisional compound chondrule formation 

We can estimate the frequency of compound chondrules 
that resulted from collisions between molten chondrules, 
assuming perfect sticking: 


/ = navct, 


where / is the compound chondrule fraction, n the number 
density of molten chondrules, a the chondrule-chondrule 
collisional cross section, Vc a characteristic collision speed 
and t the time window over which collisions lead to ccs 
fusing. Note that to avoid coalescence we require t <C r. 

For a = 0.03 cm chondrules. 


cr = 47r (2 X a) ~ 0.045 cm^. 


(14) 


At i? = 2.5 AU ( Desch[ 2007 1 , the Desch Minimum Mass 
Solar Nebula model provides a high end estimate for the 
gas density model of the Solar Nebula of pg = 2.3 x 
10“®gcm“^ at the midplane. An equal mass density in 
a = 0.03 cm, p = 3gcm“^ chondrules implies 


n ~ 6.8 X 10 ® cm 


-3 


(15) 


Chondrule-sized dust grains are difficult to concentrate 
(see iJacquet et al.||2012| iHubbard & Ebel||2015l for dis¬ 


cussions of dust concentration focusing on chondrules), so 
Equation (15) is a safe, significant over-estimate for the 


chondrule number density. 

Using disaggregated chondrules, Gooding & Keil (19811 
estimated that ~ 4% of chondrules are compound. TWai^ 


son et al. (19951 found a lower rate, ~ 2.4% by studying 


thin sections. Subsequently Ciesla et al. (20041 improved 


the 2D to 3D correction factors, finding instead 5%. We 
split the difference and use / = 0.045. However, this fre¬ 
quency is not the entire story. As noted by Gooding fc Keil| 
(|1981|) and especially by I Wasson et al.l (11995 1, compound 


chondrule do not appear to have been made by randomly 
picking two candidates from the pool of all chondrules 
and sticking them together. Radial pyroxene and barred 
olivine chondrules are heavily over-represented: combined 
they make up 97 of the 144 non-enveloping chondrules 


listed in Wasson et al. (1995)’s Table 7, despite making 
up only 10 — 13% of chondrules over all ( Gooding fc KeiH 
1981). If this result holds, it implies that those chondrule 


classes have compound chondrule frequencies higher than 
the average by a factor of over 5, and / = 0.045 is an 
underestimate of the required cc frequency. 

For viscosities below 77 < ptc (i-e. temperatures above 
T > 1095 K), chondrules collide as liquids, and can stick 


even for elevated collisional velocities. However, at T = 
1095 K, T ~ 75 s, while at only slightly hotter tempera¬ 
tures, T ~ 1145 K, T = 10 s. We require t t to avoid 
significant coalescence of the component chondrules into a 
single whole, and use t = 10 s as an overestimate, because 
it can only apply to chondrules which collide just above 
T = 1095 K. Chondrule precursors, loosely stuck together 
clumps of dust, c ould not survive coll ision speeds above 
Vc = lOOcms”^ 


(Guttler et al. 


2010 ), constraining the 
background chondrule-chondrule relative velocity. How¬ 
ever, several kms“^ shocks, one of the proposed chondrule 


formation processes (Hood & Horanyi 

1991 

Hood 

1998), 

can lead to significantly higher collision speeds ( 

Nakamoto 


(13) & Miura , 2004 Ciesla 2006). Shocks capable of reaching 


the temperatures required for chondrule formation how¬ 
ever also keep the chondrules above 1145 K for hundreds 


of seconds to tens of minutes or even longer (Holey et al. 


2013 Stammler & Dullemond, 2014). This is many times 


the chondrule stopping time in the gas of under a minute, 
so in a shock-melting scenario, once the temperature had 
fallen below 1150 K, the chondrule relative velocity would 
have fallen far below the shock’s intrinsic velocity of several 
kms“^. We therefore use Vc = 10^ cms“^ as an upper limit 
to chondrule relative velocites; while noting that restrict¬ 
ing ourselves to temperatures in the range of 1095 — 1145 K 
and time scales of t = 10 s ignores the far longer time scales 
spent at higher temperatures, which would have yielded 
significantly more collisions. Those collisions would often 


have been destructive (Jacquet & Thompson, 2014), and 


even when they led to sticking, would have resulted in co¬ 
alescence, rather than cc fusing. 

Solving Equation (13) for /, using t = 10 s, 'Cc = 
10 ^ cms“^ and the above values for n and a, we find 


/ = 0.03. 


(16) 


While this is close to the observed value, it requires a 
signihcant overestimate in the chondrule number density. 
Further, it requires that all the collisions occur at tem¬ 
peratures above T = 1095 K but significantly below T = 
1145 K while cooling out of that window on times neither 
much longer nor shorter than 10 s. Many chondrules expe¬ 


rienced multiple heating events (Jones & Danielson 1997 
Wasson & Rubin 2003 Jones et al., 2005| ), but any event 
which remelted a chondrule (T > 1700 K) would have also 
led to coalescence, leaving only the ccs that fused follow¬ 
ing the final heating event. Accordingly, we can safely dis¬ 
count molten-chondrule collisions as a cc formation mech¬ 
anism. 

For viscosities above 77 > phc (i.e. temperatures below 
T < 1095 K), chondrules collide as solids, and will only 


stick for collision velocities (Guttler et al. 2010) 


5, 0.1 cms 


(17) 


Using Equation (13) to match / = 0.045 then requires 

t~1.5xl0®s, (18) 
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or over two weeks. The requirement that t t then im¬ 
plies T < 950 K. Note further that we have overestimated 
n, and that the overall disk temperature at i? = 2.5 AU 


was much colder than 950 K (Hayashi 1981 Desch 20071. 


In the solid-collision scenario it is far easier for the compo¬ 
nent chondrules to have collided at temperatures far be¬ 
low 950 K (with almost arbitrarily long time scales avail¬ 
able) and been subsequently heated to temperatures above 
T > 900 K for just long enough for the first hints on coa¬ 
lescence to have occurred. 

5. Conclusions 

Figure shows that for any reasonably warm temper¬ 
ature, there is a fusing time scale for chondrules in con¬ 
tact to form non-enveloping ccs. If the potential cc re¬ 
mained at that temperature for too short a time, then it 
would not have fused, while if it remained at that tem¬ 
perature too long the individual components would have 
coalesced. The molten-collision model for cc formation 
makes testable predictions about the temperature of the 
chondrules upon collision (they have survived the collision 
and stuck together) and the time scales available (suffi¬ 
cient CCS must have been produced). The shaded band in 
Figure]^ shows that if “molten” chondrules are sufficiently 
inviscid that they can behave as liquids upon collision and 
stick (77 < r]bc), they are far too inviscid not to coalesce 
too fast for a sufficient number of compound chondrules 
to have been produced. This constraint does not depend 
on the precise temperature dependence of the viscosity ry, 
only the critical viscosity above which chondrules col¬ 
lide as solid, not liquid, spheres, and hence is a robust 
result. 

If compound chondrules could not have resulted from 
the collisions of molten chondrules, and indeed the pri¬ 
mary and secondary chondrules could not ever have been 
in contact at temperatures significantly above ~ 1000 K, 
then they must be the result of chondrules colliding and 
sticking while cold. While naked chondrules are not ex¬ 
pected to have been particularly sticky, experiments have 
shown that dusty rims could have bypassed that limitation 


(Beitz et al. 2012). Even so, sticking would have required 
low relative velocities, and hence the collision rate would 
have been very low, compensated for by almost arbitrarily 
long time scales chondrules could have spent floating freely 
and cold in the Solar Nebula. 

In Figure]^ we have marked the temperatures for which 
the coalescence time r is ten seconds, one hour, day and 
week, as well as the temperature below which bouncing is 
expected. The evidence that chondrules saw multiple heat¬ 


leaving a cc. The temporal longevity and strength of tem¬ 
perature fluctuations in protoplanetary disks are not well 
understood, but it seems unlikely that a proto-compound 
chondrule could be held at temperatures above ^ 1025 K 
only for time scales short enough not to coalesce (no more 
than an hour). Similarly, it seems unlikely that a proto- 
cc could be held at temperatures just above, but only just 
above, 900 K for any significant fraction of the correspond¬ 
ing T = 1 decade, so at colder temperatures there would be 
inadequate spreading to allow ccs to fuse. We also know 
from noble gas measurements that only a modest fraction 
of the solids that were eventually incorporated into chon¬ 
drites ever saw temperatures above ^ 800K (Hubbard & 


Ebel 20151, further constraining the ability to store proto- 


ccs in regions with T > 800 K for prolonged periods. We 
conclude that compound chondrules fused in a tempera¬ 
ture window of 900 — 1025 K, with very finite dwell times. 

The different viscosities expected for chondrules of av¬ 
erage composition and BO chondrules (due to the later’s 
smaller Si 02 fraction) could help explain one of the mys¬ 
teries associated with ccs: the fraction of radial pyroxene 
and barred olivine chondrules in ccs appears to be sev¬ 
eral times higher than in the overall chondrule population 


(Wasson et al. 1995). If that result holds, the participant 


pairs in ccs were not selected at random. BO chondrules 
are expected to have been less viscous than an average 
chondrule (Eigurej^, which means that they would have 
been more prone to spreading. That would have allowed 
BO chondrules to partially fuse and form ccs as a result of 
temperature fluctuations too weak to affect more viscous 
chondrules. 

The frequency of compound chondrules cannot be used 
as a probe of the chondrule density in melting regions, as 
it has been used. However, with better statistics ccs could 
potentially be used as a probe of cold chondrule-chondrule 
sticking rates. This is important as chondrule agglomer¬ 
ations are an significant part of chondrite formation, and 
there are indications that chondrules did not dwell for long 


in chondrule forming regions (Hubbard & Ebel 2015) af¬ 


ter forming. Similarly, this means that ccs provide a po¬ 
tential probe for a new temperature window in the Solar 


nebula: noble gases in trapped in chondrite matrix (Huss 


& Lewis, 1994) and primary troilite in chondrules (Rubin 


et al. 1999) place upper experienced-temperature limits of 


about 800 and 650 K respectively, significantly colder than 
our 900 — 1025 K. This new temperature window is espe¬ 
cially interesting as it lies near the ^ 1000 K threshold for 
thermal ionization to allow the magneto-rotational insta¬ 


bility to operate (Gammie 1996). These temperatures can 


also be further constrained by considering isotopic frac¬ 


ing events (Jones et al. 2005) implies that many such pairs tionation or lack thereof ( Alexander et al. 2000) because 


of stuck together (but not yet fused together) chondrules 
would have experienced episodes of elevated temperature. 
If those episodes were a small but still respectable fraction 
of the coalescence time of their peak temperature, the com¬ 
ponent chondrules would have fused, but not coalesced. 


our temperature window is close to the 50% condensation 
temperature of potassium in the Solar Nebula (^ 1006 K, 


Lodders 2003). 
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